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RATIONAL TRIGONOMETRIC APPROXIMATIONS 
USING FOURIER SERIES PARTIAL SUMS 


James F. Geer 1 

Department of Systems Science 
Watson School of Engineering and Applied Science 
State University of New York 
Binghamton, NY 13902 


ABSTRACT 

A class of approximations { *Sat,a/ } to a periodic function / which uses the ideas of Pade, or 
rational function, approximations based on the Fourier series representation of /, rather than 
on the Taylor series representation of /, is introduced and studied. Each approximation Sj\\m 
is the quotient of a trigonometric polynomial of degree N and a trigonometric polynomial 
of degree M. The coefficients in these polynomials are determined by requiring that an 
appropriate number of the Fourier coefficients of N/v,m a S ree with those of /. Explicit 
expressions are derived for these coefficients in terms of the Fourier coefficients of /. It is 
proven that these “Fourier-Pade” approximations converge point-wise to (f(x+) + f(x ~))/ 2 
more rapidly (in some cases by a factor of \/k 2M ) than the Fourier series partial sums 
on which they are based. The approximations are illustrated by several examples and an 
application to the solution of an initial, boundary value problem for the simple heat equation 
is presented. 


! This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NASI- 19480 while the author was in residence at the Institute for Computer Applications in 
Science and Engineering (IOASE), NASA Langley Research Center, Hampton, VA 23681. 
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1. Introduction. Fourier series are used widely in many branches of applied math- 
ematics. For example, they are often used together with separation of variables to con- 
struct analytical solutions to boundary value problems for differential equations and 
with a variety of spectral methods to find approximate solutions to these problems nu- 
merically. For practical purposes, approximate solutions to these problems are often 
obtained using only a finite number of the terms in a Fourier series. This trunca- 
tion procedure often leads to nonuniformly valid approximations. For example, when 
the function being approximated has a point of discontinuity, the Gibbs phenomena is 
present. The “oscillations” caused by this phenomena typically propagate into regions 
away from the singularity, and, hence, degrade the quality of the partial sum approx- 
imation in these regions. Even if the function being modeled is analytic, but has a 
region of large slope or curvature, there can be significant oscillations in the partial 
sums outside these regions, which again degrade the quality of the approximation. 

Recently, Gottlieb and Shu [6] and Gottlieb, et.al. [7] have proposed a way of 
overcoming the Gibbs phenomena. Their technique involves the construction of a new 
series using the Gegenbauer polynomials C A (x). For a function / that is analytic on 
the interval [-1,1], but is not periodic, they prove that their technique leads to a series 
which converges exponentially to / in the maximum norm. To do this, they require 
that the parameter A, which appears in the weight factor (1 - x 2 3 ) A_1/2 , grows with 
the number of Fourier modes considered. As we shall demonstrate below, the family 
of approximations we shall introduce can be defined explicitly in terms of the known 
Fourier coefficients. This will prove to be particularly useful for certain applications (as 
we shall demonstrate) when the Fourier coefficients are themselves functions of one or 
more other variables. Although the approximations we shall define do not “eliminate” 
the Gibbs phenomena, they do mitigate its effect, as we shall show. This is especially 
true outside a “small” neighborhood of a point of discontinuity of /, where, for practical 
purposes, the “unwanted” oscillations can essentially be eliminated. 

To fix notation, we let / (a:) be a piece-wise smooth, 27r-periodic function. Then 
we can associate with f(x ) its Fourier series S(x) defined by 

OO 

(1) S(x) = a 0 / 2 + ^(a n cos(nx) + 6„sin(nx)), 

n=l 

(2) a n = ( 1 / tt ) /: f(x)cos(nx)dx, b n = (l/?r) J f(x)sin(nx)dx, 

n — 0 , 1,2, ... . (For some of our formulas below, it is convenient to regard the coefficients 
{a n , b n } as being defined by (2) for negative as well as positive values of n. Thus, a_ n = 
a n and 6_ n = -b n for all integers n.) It is well known (see [2], for example) that S(x) 
converges to f(x) at each point x where / is continuous and to (/(x + )-f f(x~))/2 at each 
point x where / is not continuous. For practical purposes, f(x) is often approximated 
by its Fourier partial sum Sn,o{x) defined by 
N 

(3) Sn,o(x ) = a 0 / 2 + E<«. cos(nx) + b n sin(nx)). 

n— 1 
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We now define a class of approximations to / which uses the ideas of Pade approx- 
imants (see, e.g., [1]), except that the approximations are based on the Fourier series 
representation of /, rather than on the Taylor series representation of /. Thus, for any 
two non-negative integers N and M, we define a family of “Fourier-Pade” approxima- 
tions SN'M(x) by 


,.s C ( \ A 0 /2 + '£ i * =l (A n cos{nx) + B n sm(nx)) 

1 4 1 vnm (x) = M n 7 r — 77 — t~t — rr - 

1 + Em=i (^m cos(mx) + D m sinfrax)) 

The 2 M + ‘2N + 1 constants {A n , B n , C m , D m } which appear in the definition of 
are determined by the condition that 2 M -f 2 N + 1 of the Fourier coefficients of Sn,m 
agree with those of /, i.e. 


(5) 


(l/*) J Sn,m( X ) cos (jx)dx = a 3 , 0 < j < J 


(6) (l/ 7 '') f Sn,m( x ) sln(kx)dx = 6^, 1 < k < K 

J —IT 

where J + K = 2 M -f 2 N. 

The idea of constructing Pade approximations based on series representations of 
functions other than the classical power series representation has been suggested and 
studied by several other investigators. For example, Maehly [9] has suggested an ap- 
proach to determine the coefficients in rational approximations based on Chebyshev 
series, and this approach has since been described in numerical terms, with examples, 
by Ralston[lO] and by Fike[5] . An account of the basic theory of such approximations 
has been given by Cheney [3], who considered the more general case of expansions in 
terms of a basis {ifj}, whose elements satisfy relations of the form <pi<?j = a ijk i Pk- 
Clenshaw and Lord [4] have reviewed rational approximations based on Chebyshev series 
and present a well-conditioned method for constructing the elements of a “Chebyshev- 
Pade” table. The general form of a Fourier-Pade representation, such as equation (4), 
has also been suggested by Cheney. He has proven the existence of a rational trigono- 
metric function which “best” approximates a continuous function /, but does not discuss 
any detailed algorithms for the computation of the coefficients, rates of convergence of 
the approximations, etc. 

In the sections 2-4 below we consider the special case of equation (4) when M = 1. 
In particular, in sections 2-3 we develop explicit formulas for the coefficients which 
appear in the definition of Sn , i in terms of the Fourier coefficients of /. Several results 
concerning the behavior and convergence of the approximations {6^,1} are presented in 
section 4. In sections 5-6, the results of sections 2-4 are generalized to values of M > 2. 
In section 7 the approximations are applied to a simple initial, boundary- value problem 
for the heat equation. Some observations and insights about this class of approximations 
are discussed in section 8. 
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2. Odd functions of x. Since / can be expressed as the sum of an even function 
of x and an odd function of x , it is sufficient to consider separately the cases when / 
is either even or odd. Thus, in this section, we consider the case when / is an odd 
function of x, while the case of an even function will be discussed in the next section. 
For this case, the Fourier series of / can be expressed as (1) with each a n — 0, n > 0, 
and we define a class of Fourier-Pade approximations Sn,m( x ) °f the form 


(7) 


Sn,m( x ) 


En=i Bn sin(nx) 

1 + Em=l C 'mCOs(mx)' 


In this section we examine in detail the case M = 1. The case M > 2 will be discussed 
in section 5. Thus we consider first approximations of the form 


(8) 


O / x E«=i B n sin(nx) 
= 1 + Ci cos(x) • 


The constants {B x , B 2 , . . ■ , Bn, Ci} are determined by conditions (6), which for this 
case become 


N 


(9) 


where 


T.dkjBj = bk, k = l,...,N+l, 

3=1 


, /, / \ f* sm{jx)sm{kx) 

^ = C /*>/_, TTc^(x) iz = ii 


To solve (9), we first use the identity 


sin(jx) sin (kx) — (l/2)(cos((A: — j)x) — cos ((k + j)x) 


to express the coefficients {d^j} as 

tit COS! TlX ) 

(10) d k „ = I lt . M - /**„ where = (1 /t) ^ Q ^dx, n > 0. 
Assuming that 0 < C 2 < 1, we find (see, e.g., [8], p. 113) 

2p 


1 + o 2 \fi ~ C\ - 1 

(11) /„,i = p n - -, for n > 0, where p = 

1 — p* C-i 


or Ci = 


1 +p 2 


Using this result in (10) we can write 

d k ,j = p k ~ j ( 1 + p 2 )(l + P 2 + + p 2l ~ 2 ), for k > b 

Using these expressions, (9) can be expressed in matrix-vector form as 


( 12 ) 


5ft 


B\ \ 


bl \ 

b 2 

1 

b 2 

\ Bn ) 

1 + p 2 

&N+1 / 
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where the (N + 1) x N matrix 3? is defined by 


1 P P 2 

P 1 + P 2 P{^+P 2 ) 

p 2 p(l + p 2 ) 1 + p 2 + p 4 


p N - 1 \ 

p N ~ 2 (l+p 2 ) 
p N ~ 3 ( 1 +p 2 + p 4 ) 


p N ~ 1 

^ P N ' 1 (^ + P 2 ) 


1 + p 2 + --- + p 2N ~ 2 
p(l + p 2 + ■ ■ ■ + p 2N ~ 2 ) ) 


To solve (12), we first note that if we subtract p times the next to last equation from 
the last equation we obtain the relation 0 = (&w+i — P^n)/{ 1 + p 2 ), from which we find 
(assuming b jv ^ 0) 


(13) 


P = 


fyv+i 

bs 


or Ci = — 


2p 


l+p 2 


26jyfrjV+i 

+ ^AT+l 


In a similar manner, subtracting p times the {k — l) </l equation from the k th equation 
we find 


Bk = ~ TXT 1 - E P k ' J B 3 , l < k < N. 

1 + P j=k + 1 

Using this recursive definition of the {2?*}, we can use induction to show that the {5*} 
are given explicitly by 

(14) B k = b k - ^hzL± b Jl±}l =bk + °l{b k _i + b k+i ) 1 <k<N. 

1 + p l 2 

(We note that equations (14) also follow from equating (8) to bj sin(jx), multiply- 

ing by 1 + Ci cos(x), and equating coefficients of sin(nx) on each side of the resulting 
expression, for 1 < n < N.) Thus, the coefficients {2?*} and C\, which appear in 
the definition (8) of Sn, i, are completely defined in terms of the coefficients {&*,} by 
equations (13) and (14). 

Example 1: As an example of these results, we let f(x) = x/2, for — tt < x < ir, 
and f(x + 2x) = f(x). Then b k = ( — 1 ) k+l /k, for k = 1,2,... . Using the formulas 
above, we find p = — N/(N + 1) and hence 

22V(7V + 1) _ 32V(2V + l) + 2 

1 2N{N + 1) + 1’ 1 4N{N + l) + 2' 


B k = (~l) k 


2N(N + 1) k 
2N{N + 1) + 1 k* - 1 


1^, = 2,3, ... ,2V. 


Explicit values for C\ and B\, ...,2?jv for a few small values of TV are shown in Table I. 
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TABLE I ( Example 1: M = 1) 


N 

Ci 

Bi 

b 2 

S3 

B, 

b 5 

Be 

1 

4/5 

4/5 






2 

12/13 

10/13 

3/26 





3 

24/25 

19/25 

7/50 

-2/75 




4 

40/41 

31/41 

37/246 

-4/123 

5/492 



5 

60/61 

46/61 

19/122 

-13/366 

3/244 

-3/610 


6 

84/85 

64/85 

27/170 

-19/510 

23/1700 

-1/170 

7/2550 


To illustrate these results, in Figure 1 we have plotted the function /(x) (solid line), 
along with the approximations £4,0(2:) (dashed line) and £3,1 (x) (dotted line). As the 
figure illustrates, except possibly in a neighborhood of the discontinuities of /, £3,1 is 
a better approximation to / than is £4,0. In particular, the oscillations present in £4,0 
away from the singularities of / have been virtually eliminated in £3,1. However, near 
the points of discontinuity of /, the Gibbs phenomenon is still present in £3,1 , although 
its effect appears to be somewhat less than the corresponding phenomena present in 
£4,0- We shall discuss this example further in sections 4 and 8. 

Before continuing, we briefly note two special cases. First, if 6 jv+i = 0, then p = 0 
and hence Ci = 0. Thus Bk — bk, for 1 < k < N, and £jv,i reduces to the Fourier partial 
sum £jv, 0, for this case. Second, suppose 6jt = 0 when k is even. In this case, we modify 
the form of £;v,i and define 


En=i B 2n -i sin((2n - l)x) 

1 + C 2 cos(2x) 

Then equations (6) for this case have the same form as (9) with bk replaced by b 2 k-i 
and dk,j = d 3 k = I\k-j\ ~ h+j-u from which we find 


(15) 


£jv,i(s) 


dkJ = for * > j. 


Using these expressions and the same type of arguments as presented above, we find 
for this case (assuming 6 2 jv- 1 7^ 0) 


(16) 


C 2 = 


_ l>2N+l 
1 + p 2 ' P b 2 N-\ ’ 


C 2 

B 2 k-i = b 2 k-i + —{b 2 k-3 


+ ^2Jt+l)i 1 £: k < N. 


(Recall that we have defined 6_i = —61.) 
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3. Even functions of x. When / is an even function of x, its Fourier series has the 
form (1) with each b n = 0. Then we define the family of Fourier- Pade approximations 


(17) W*) = 

1 + £m= 1 C m COs(mx) 

In this section we consider the case M = 1 . The case M > 2 will be discussed in section 
5. Thus we consider approximations of the form 


Sn,i{x) = 


A/ 2 + £n=l An cos(nx) 
1 + C\ cos(x) 


The requirement that the first N + 2 Fourier coefficients of Sn i \ agree with those of / 
leads to a system of equations for the TV + 2 unknowns {Ci, A 0 , ..., Ayv), which can be 
solved using the techniques of section 2. In particular, we find (assuming a n ^ 0) 

(19) Ci = - y , where p = — A k = a k + ^-(afc-i + a k+l ), 0 < k < N. 

1 + p 2 a N 2 

In the special case when ajv+i = 0, we see that p = 0 and hence Ci = 0. Then each 
A k = a k and i(x) = 5jv,o(®) for this case. Also, in the special case that a k = 0 when 
k is odd , we redefine 


Sn,i(x) = 


A )/ 2 + En=i A 2n cos(2wx) 
1 + C 2 cos(2x) 


and find (assuming a 2 v 7 ^ 0 ) 


(21)C 2 = , where p - 2^+1 ^ ^ - a2k + ^(a 2fc _ 2 + a 2Jt+2 ), 0 < k < N. 

It p CI2N & 

Example 2: As an application of these results, we let /(x) = 1/y^l + acos(x), 
where a 2 < 1. Here / is analytic for — tt < x < ir, but develops a sharp, narrow peak 
near 1 = 0 as a - >— 1. In Table II we have recorded the Fourier coefficients {a n } of /, 
as well as the coefficients {Ci, Aq, ..., Ajv} for a few small values of N when a = —0.95. 


Table II ( Example 2: a = —0.95, M = 1) 


N 

on 

c, 

Aq 

A\ 

A 2 

Aq 

a 4 

0 

2.94734 

-0.68380 

2.15059 





1 

1.16518 

-0.85289 

1.95357 

-0.37012 




2 

0.65289 

-0.89142 

1.90867 

-0.43948 

-0.04495 



3 

0.40050 

-0.90812 

1.88922 

-0.46953 

-0.05802 

-0.01235 


4 

0.25636 

-0.91741 

1.87839 

-0.48626 

-0.06529 

-0.01657 

-0.00452 
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To illustrate these results, in Figure 2 we have plotted the function f(x) (solid 
line), along with the approximations S4 i o(x) (dashed line) and 53 ,i(x) (dotted line). As 
the figure illustrates, .S^i is a noticeably better approximation to / than is ^4,0. In 
particular, the oscillations present in S^o have been virtually eliminated in . We 
shall discuss this example further in sections 4 and 8. 

4. Analysis of the Case M = 1. We now prove certain results which express 
the asymptotic behavior of the coefficients C\ and {A*} or {Bk} as N — > 00, and which 
also allow us to make some statements concerning the convergence of the family of 
approximations {S^i}, as N — > 00. 

It is well known (see [2], for example) that if / has q continuous derivatives on 
— 7r < x < 7r and the derivative of / of order q + 1 is piece-wise continuous, then the 
Fourier coefficients of / are 0(1/ k q+2 ) as k — ► 00. For such a function, we now show 
that the coefficients in Sn, 1 decay as k — ► 00 at a faster rate than the Fourier coefficients 
on which they are based. 

THEOREM 1. Let f be an odd, 2n-periodic, piece-wise smooth function and let its 
Fourier coefficients {&*} satisfy the condition that bk = 0(l/k p ), as k —* 00, where p is 
a real, positive number. Then \C\\ = 1 + 0(1/N 2 ) and Bk = 0(1/ k p+2 ) as k,N — * 00. 
More precisely, let b, (3, and~f be constants (independent of k). Then, ifbk — (b/k p )( 1 + 
jd/k - |- 7 jk 2 + 0(1 jk z )), as k — ► 00, it follows that 

C ' = + VO + 0(1/A,3) ' Bt = J i (is;) +p(1 - (t/JV)!) ) + 0( ST5>* 

as k,N —y 00. If f is an even, 2n-periodic, piece- wise smooth function, then these results 
hold with bk replaced by ak and Bk replaced by Ak . 

Proof: We shall outline the proof only for the case when / is an odd, 27r-periodic, 
piece-wise smooth function, since the proof for the case when / is an even function 
follows the same line of reasoning. Using the assumed form of the asymptotic behavior 
of the coefficients {&*}, we note first that we can write 

V + . = (V(JV + l)')(l + /3/{JV + l) + 7 /(JV + l) ! + 0(l/(7V + lf)) 

(22) = (l/iV)(l+£)-f(l+|9 e/(l+e) + 7fy(l+e) 2 + 0(l/Ar 3 )) 

= (b/W )(1 + (0 - p)e + (7 - fi(p + 1) + p(p + l)/2) £ 2 + Ofl/JV 1 )), 

where we have defined e = 1/N. Then, using this expression and the definition of p in 
(13), we can write 

< 23 > ' = ^ = 1 -^ + e± ^ M + °( , / jv3 )'“ w -~- 
and hence 

Cl = -7^5 = -1 + ff, + 0(l/N% 

as stated in the Theorem. 
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In a similar manner, for large values of k we can use the assumed form of the 
coefficients {bk} to write 

b k ± i = (b/(k±l) p )(l+/3/(k±l) + 1 /(k±l) 2 + 0(l/F 3 )) 

= (b/k p )( 1 ± e)" p (l + 0t/( 1 ±e) + 'ye 2 / ( 1 ± e) 2 + 0(1/ k 3 )) 

= (ty*-)(l + (/? T P)e + (7 =F «P+ 1) +p(p + l)/2)e ! + 0(1/* 3 )), 

where we have now defined e = 1/A:. Using these expressions, along with an analogous 
expression for bk and the asymptotic expression above for C\, in the definition of Bk in 
(14), we find for large values of k and N (with k < N) 

Bk = bk 4- (C\/2)(bk-i + frjt+i) 

= (6/A; p ){l + /3e + ye 2 + 0(e 3 ) 

+(l/2)(— 1 + p 2 /(2N 2 ) + 0(1/N 3 ))(2 + 2ffi + (2 7 + p 4- p 2 )e 2 + 0(e 3 ))} 

= —(bp/2)(l/k p+2 )(\ +p( 1 - (k/N) 2 )) + 0(1/ k p+3 ), 


which completes our proof. 

Thus, Theorem 1 establishes that the coefficients {Z? n } in the approximations Sjv.i 
decay to zero more rapidly (by a factor of k~ 2 ) than the coefficients in the partial sum 
Sn , o on which they are based. 

Before continuing, we remark that results very similar to those stated in Theorem 
1 hold if the coefficients {6*} of / have the same asymptotic form as that stated in the 
theorem, except for a multiplicative factor of (— l) fc , i.e., if bk = (b/k p )( — l) fc (l + fl/k + 
'y/k 2 + 0(l/k 3 )), as k — ► oo. In this case, it is easy to show that the Fourier coefficients 
{6jt} of / = f(x + 7r) are related to those of f by bk = ( — l) k bk- It then follows that 
the coefficients {bk} satisfy the conditions of the theorem and hence the corresponding 
coefficients C\ and { B k } have the asymptotic form indicated in the theorem. In partic- 
ular, it follows from equations (13) and (14) that C\ = — C\ and Bk = (— 1 ) k B k - Thus, 
the asymptotic behavior of C\ is the negative of that indicated for C\ in the theorem, 
while the asymptotic form of B k is (— l) fc times the form indicated for B k - 

To illustrate these results, in Figure 3 we have used the coefficients {Bk} from 
Example 1 and have plotted fc 3 |Bfc| as a function of 1/k, for 3 < N < 20. The figure 
clearly illustrates the 0(l/k 3 ) decay of the coefficients {Bk}, and also indicates an 
interesting asymptotic behavior of B^j- q , for a fixed value of q, as N —* oo. 

If the function / is analytic for — n < x < n, then (see, e.g., [2]) the Fourier 
coefficients of / decay exponentially as k — > oo. That is, there exists a constant 9, with 
0 < |0| < 1, such that are 0(9 k ) as k — ► oo. The next theorem shows that, 

again for this case, the coefficients {Bk} decay to zero more rapidly than the Fourier 
coefficients on which they are based. 

THEOREM 2. Let f be an odd, analytic, 2n -periodic function and let its Fourier 
coefficients { 6 *,} satisfy the condition that b k = 0(\8\ k / k p ), as k — » oo, where 0 < \8\ < 1 , 
and p is a real number. Then C\ = 29/(1 + 9 2 ) + 0(1 /N) and Bk = 0(\9\ k / k p+l ) as 
k,N — ► oo. More precisely, let b, /3, and 7 be constants (independent of k). Then, if 
bk = (b9 k /k p )( 1 + fl/k + 'y/k 2 + 0(1 /k 3 )), as k — > 00 , it follows that 


1-P- 

N 


1 -9 2 
1 +8 2 


2/?)( 1 — 9 4 ) + p 2 (l — 60 2 + 8 4 ) 
2(1 + 8 2 ) 2 N 2 


+ 0(1/N 3 ) 
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b„ = U^d - */*) - j 5rr^ (2 ^ + (W + p)(1 + k/N))(1 ~ k,N) 


P 

2(1 + 9 2 ) 2 


(d + « ! ) 2 - am i - 6« 2 + o) (i - */jv) - } + o(-^j), 


as k,N —* oo. If f is an even, analytic, 2n- periodic function, then these results hold 
with bk replaced by ak and Bk replaced by Ak- 

Proof: As in the proof of Theorem 1, we shall outline the proof only for the case 
when / is an odd function, since the proof for the case when / is an even function 
follows the same line of reasoning. Using the assumed form of the asymptotic behavior 
of the coefficients {£»*,}, we note first that b^+i has the asymptotic form indicated on 
the right side of equation (22), multiplied by 9 N+1 . Then, using this expression in the 
definition of p in (13), we find that p has the asymptotic form indicated by the right 
side of (23) multiplied by 9. Using these expressions in the definition of C\ in (13) we 
find 


Ci = 


1 +p 2 


= -29 


1 - P/N ± (p ± p 2 - 2/?)/(2iV 2 ) + Q(l/iV 3 ) 

1 + 9* (1 - p/N + {p + p 2 - 2f3)/(2N 2 ) + 0(1/N 3 )f 


29 j p /l-fl 2 \ . (p ~ 2/3)( 1 — fl 4 ) + p 2 (l - 6 fl 2 + (P) 0(1/m \ 

~ 1+0 2 ! 1 N \ 1 + 0 2 J + 2(1 + 9 2 ) 2 N 2 + U j J’ 

as AT — ► oo. 

In a similar manner, using the assumed asymptotic form of the coefficients { 6 ^} we 
can write 

b k ± i = (W fc± 1 /F)( 1 + (^p)e + ( 7 T^+l) + p(p + l)/ 2 )e 2 + 0(1/ k 3 )), 
as k — > oo, where e = 1/k. Thus, for large values of both k and N we have 

Bk = bk + (C\/2)(b k -\ + bk+i) 


2 1 fi P . (p-2/9)(l-0 4 ) + p 2 (l- 6 ^ 2 +^) 

7e l+0 2 [ N\l+9y + 2(1 + 9 2 ) 2 N 2 


■ 1 + (/3 + p)e + (7 + f}(p + 1 ) + p(p + l)/ 2 )e 2 
+ & (l + {P ~ p)e + (7 - P(p + 1) + P(P + l)/2)e 2 )] } + 0(9 k /k p+3 ) 

= ~£i {pj^d - W - i [5^(2^ + w + P)d + W>d - k/N) 
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-2fIT vy (<> + ** - (*/W - ^ + *)) (. - k/N) - } + 0(^3), 

as k — * oo. This completes our proof. 

We note that the rate at which the coefficients {£?*} decay to zero changes slightly 
as k —> N. In particular, for k < N, Bk = 0(O k ( 1 — 0 2 )(1 — k/N)/k p+l , but B n = 
O(0 N /N p+2 ), as k, N — > oo. Thus, while each Bk decays to zero faster than the corre- 
sponding Fourier coefficient bk, we see that Bk , decays somewhat faster than Bk 2 , when 
k 2 < ki < N. 

To illustrate these results, in Figure 4 we have used the coefficients {A*} from 
Example 2 and have plotted A: 2 |y4* |/|a*| as a function of 1 jk, for 3 < N < 20. The 
figure clearly illustrates that the coefficients { Ak } decay at a faster rate (asymptotically 
by a factor of k~ 2 ) than the coefficients {a*}, and also indicates, as in Example 1, an 
interesting asymptotic behavior of Ajv_ 9 , for a fixed value of q, as TV — ► oo. 

Finally, we use the results of the previous two theorems to show the manner in 
which the approximations {Sjv,i(#)} converge to the original function f(x). 

THEOREM 3. Let f be an odd, 2 tt - periodic, piece-wise smooth function and let 
its Fourier coefficients {&*} satisfy the conditions of Theorem 1. Then the sequence of 
approximations (5jv,i(x)} converges to ( f(x + ) + f(x~))/ 2, as N — > oo, for all —w < 
x < 7r. Moreover, the sequence {5jv,i(*)} converges like a Fourier series whose terms 
are 0(l/iV p+2 ). These results also hold if f is an even, 2n -periodic, piece-wise smooth 
function, if the Fourier coefficients {a*,} of f satisfy the conditions stated above with bk 
replaced by ak . 

If f is an odd, analytic, 2ir -periodic function and its Fourier coefficients satisfy the 
conditions of Theorem 2, then the sequence of approximations {5^,1 (a;)} converges to 
f(x), as N — > oo, for all —it < x < n. Moreover, the sequence {5jv,i(a0} converges like 
a Fourier series whose terms are 0(9 N /N p+1 ). These results also hold if f is an even, 
analytic, 2ir -periodic function, if the Fourier coefficients {ak} of f satisfy the conditions 
stated above with bk replaced by ak. 

Proof. As in the proofs of the previous theorems, we shall outline only the proof 
for the case when / is an odd, 27r-periodic, piece-wise smooth function. 

To begin, we note that, since / is an odd function, we can restrict our attention to 
values of x in the open interval 0 < x < it (since Sjv,i(0) = (/(0 + ) + /( 0 — ))/2 = 0 = 
Sn,i(k) = (/(^ + ) + /( 7r ~))/ 2, and hence convergence is assured at x = 0 and x = ir in 
a trivial way). Consequently, we let x be any fixed number in the interval 0 < x < rr 
and let c > 0 be any fixed positive number. We then define 

(24) £*.,(*) = (/(s + ) + /(*-))/ 2 - S N .,(x) = £$>(*) + £#>(*), 

(25) Eg\z) = (/(*+) + /(*-))/: 2 - S w ( i), E%\x) = S N , 0 (x) - S w (i). 

By the hypotheses of the theorem concerning /, the Fourier partial sums SN t o(x) 
converge to ( f(x + ) + f[x~))/ 2, as N — * oo. Hence, there exists a positive integer 
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Ni — Ni(x,e) such that 


(26) \EW(x)\ = | (/(x + ) + /(*-))/ 2 - ^, 0 (x)| < -, for all N > N t . 

Also, using the definitions of Sn,o and Sn,u we can write 


<(*) = S K , o(x) - £„.,(*) = E 4„ sin(nx) - 


(27) 


Ci 6 at sin((iV + l)x) — 6 jv+i sin(Nx) 
2 1 + Ci cos(x) 


Using the assumed asymptotic form of the coefficients {&„}, there exists a positive 
integer N 2 = N 2 (x, e) > (4 |6| /(e(l — cos(x))) 1/,p such that \b n \ < 2 |6| /N p and C\ > — 1 
for all N > N 2 • Then, using equation (27), we can write 


\e^(x)\ = \S n , 0 (x)-S na (x)\< 


|Ci| |6at| + |6yv + i 

2 1 — cos(x) 


(28) 


< 


4 1 *| 


< for all N > N 2 . 


2(1 — cos(z))AT p 2’ 

Using the bounds (26) and (28), from equation (25) we find 

(29) |(/(x+) + /(x-))/2 - S Wl ,(x)| < |Bg’(x)| + |£{?<x)| < | + j = £, 

for all N > N 3 = Max(N\, jV 2 ). Thus, the sequence of approximations { 5/v,i } converges 
to (/(x + ) + /(x~))/2 as N — ♦ 00 . 

To demonstrate the manner in which these approximations converge, we write Sjy t % 
as 

N 

(30) S N , 1 = where r x = 5i,i and r k = S k ,i(x) - S k - i,i(x), 2 < k < N, 

k= 1 

and then examine the quantities r k for large values of k. Using the definition (8), we 
can write 


(31) 


r k (x) = 


En=i BM sin H En=i sin(n:r) 

1 + Ci W cos(x) 1 + cf fc_1 ^ cos(x) 


where we have placed a superscript on the coefficients B n and Ci , to remind us of the 
dependence of these quantities on the parameter k. Combining the terms on the right 
side of equation (31) and using various trigonometric identities, we can write 


rjt(x) = 


Num 
Den ’ 


(32) 
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where 


and 

( 33 ) 


Num = 


sin((A; — 1) 


^-C\ k) | 


bk 4 (bk-i + &jt+i) 


) 


!i ^( 6 ‘ + ^F- 6 ‘-' + T- 6 ‘«) 


I 

sin((t + l)i) c , t -i) { fct + (it - + tw) [ , 


Den = (1 + C[ k) COS (*))(! + Cf‘-" cos(x)). 


Using the asymptotic form of the coefficients B and C [ k * from Theorem 1, we find 


Num = — 


bp 1 

lfib+2 


(1 — cos(x)) sin(£x) + 0 (l/k p+3 ), 


v 2 

( 34 ) Den = (1 — cos(x)) 2 + 7 ^ (1 — cos(x)) cos(x) + 0 (l/k 4 ), as k —* 00. 

AC* 

Using the estimates ( 34 ) in ( 32 ), we find 


( 35 ) 


rk 


bp 1 sin(fcx) 
2 1 — cos(x) k ? + 2 


+ 0 (l/k p+3 ), 


as k 


oo. 


This completes the proof of our theorem. 

To illustrate some of these results, in Figure 5 we have plotted log |.E;v,i|/ log |iV| 
for Example 1, as a function of 1/ log \N\, for 3 < N < 50 , with x = 7r/2. We note that, 
if En , i ~ aN~ p , as N — + oo, then log |£;v,i|/ log \N\ ~ — p + log |a|/ log \N\, as N — ► oo. 
Thus, the plot should be approximately linear for small values of 1/ log \N\, with an 
intercept of — p. The value of p = 3 , as predicted by Theorem 3, is clearly suggested by 
the figure. In this figure we have also plotted the corresponding error for the Fourier 
partial sums, corresponding to M = 0 . For M = 0, the figure suggests a value of p = 1. 
Thus, the improvement in the rate of convergence suggested by the figure is consistent 
with the increase predicted by Theorem 3 , i.e., an increase in p of 2. 


5. Case M > 2. We now generalize the results of the previous sections. In par- 
ticular, when / is an odd function of x, we consider the family of approximations 


Sn,m( x ) 


E«=i ffnsin(nx) 

1 + Em=l C m COS (mi)’ 
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( 36 ) 



where M is any positive integer. The constants { B n ,C m } are determined by the con- 
ditions (6), which for this case become 


^2 dkjB-j = bk, k=l,...,N + M, 

3 = 1 


where 


(38) d kj = (1/x) - 


n s'm(j x) s'm(kx) 


+ £m= 1 C m COS (mi) 


dx = dj t k = I\k-j\,M ~ h+ h M, 


In,M = (1/*) jf J 


cos(nx) 

+ £m=l Cm COs(mi) 


dx, n > 0. 


To evaluate /„,m, we denote the denominator of the integrand by P and express cos(mi) 
as a polynomial of degree m in cos(x). Then P becomes a polynomial of degree M in 
cos(x) and, assuming Cm ^ 0, we can write 

M M 

(40) P = 1 + ^2 Cm cos (mi) = 2 m ~ 1 Cm jQ(cos(x) + z , ), 

m= 1 t= 1 

and hence 

1 i M a, M 

^ 41) 1 T £m=i C m cos(mx) _ 2 m ~ 1 Cm cos(x) + Zi' a * Zt) 

Here { — Z{] are the roots of the polynomial P, when P is regarded as a polynomial in the 
variable z = cos(x), and we have assumed that the { 2 ,} are all distinct. Substituting 

(41) into (39) and using (11), we can write 


9 M 

Z a Pi 


In } M = ~~p 51^*1 * "Zfl 

Cm 1 ~ Pi 


a= n 


k=\ k *i (/»* - - p«-p*) 


, pi = \JA - 1 - *»• 


Using this result in (38) we find 


0 m j - 1 


t=l /=0 

To solve equations (37), we first define the quantities 

s 0 = l, s m = (-l) m £ \<m<M. 

1<»1 <12 
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(44) 



(We note that s TO is just the coefficient of p M ~ m in the polynomial n^it/ 7 — P > )•) We 
then denote the k th equation in (37) by eq(k) and form the linear combinations 

M 

(45) lc(k) = X s m eq(k - m), k = N + M, N + M - 1, 1. 

m=0 

Here we define e<?(0) to be the trivial equation 0 = 0, and eq(n) — —eq{—n), for n < 0. 
For k = N + M, N + M — 1 , N + 1, we find that the left side of each lc(k) vanishes. 
This leads to the relations 

M 

(46) X s mbN+M+l-p-m = “ bpf+M+l-p, P = 

771=1 

which are a system of M linear equations for the M quantities We shall 

assume that the determinant of the matrix of coefficients in (46) is nonzero so that 
these equations can be solved uniquely for Once these quantities have been 

determined, we can use equations (44) and (42) to express each s m in terms of the 
quantities z,, and then use equation (41) to express each of the coefficients C m in terms 
of the z t . We write out these equations for a few small values of M and then use induction 
to show that the coefficients {C m } can be expressed explicitly in terms of the known 
quantities {s t } as 

c\ M —m M 

(47) c m = — X s kS m +k, 1 < m < M , D = X 5 1 

U *=0 j = 0 

To determine the coefficients for k = N y N — 1,..., 1, we solve lc(k) for Bk 

and then use induction to show that the coefficients {Bk} are given explicitly by 

i M 

(48) B k = b k + - X C m (£>jt +m + bk- m ), 1 < k < N. 

* m—1 

(Equations (48) also follow from equating (36) to bj sm(j x) , multiplying by 1 + 

J2m~i Cm cos(mx), and equating coefficients of sin (nx) on each side of the resulting 
expression, for 1 < n < N.) Thus, the coefficients { C m } and {B n } are completely 
determined by equations (48), (47), and (46). 

Example 3: As an application of these results, we consider the odd, 2k— periodic 
function / defined by f(x) = tx(n — x)/(e 2 + (x — S) 2 (x — n + 8 ) 2 ) ! / 2 for 0 < x < ir, 
where e and 8 are real, positive parameters. For this function, we find bk = 0 when k is 
even and hence = X)£=i &2fc-i sin((2 k — l)x). Then, for M = 2, we define 

B 2 k-i sin((2 k - l)x) 

c jfc==l 

' V ' 2 1 + Ci cos(2x) + C 2 cos(4x) ’ 

where the coefficients {H2fc-i} and {C TO } are defined by equations (46)-(48), with bj 
replaced by 62^-1. Explicit values for {Ci, C2, B\, ..., H2/V-1} are shown for a few small 
values of N in Table III when t = 0.1 and 8 — 0.5. 
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Table III ( M = 2) 


N 

b 2N-1 

Ci 

c 2 

B x 

B 3 

B 5 

b 7 

b 9 

1 

0.41814 

-.99477 

0.29067 

0.43574 





2 

0.31809 

-1.3981 

0.46568 

0.43655 

-.16521 




3 

0.09684 

-1.3559 

0.54610 

0.42554 

-.17564 

0.02412 



4 

-.11150 

-1.3462 

0.58487 

0.42077 

-.18343 

0.02989 

0.00716 


5 

-.17153 

-1.3539 

0.59107 

0.42047 

-.18703 

0.02987 

0.00817 

0.00234 


In Figure 6, we have plotted 5g,o» 5s, l (using equation (15)), and 54,2, as well as /, for 
0 < x < 7T. From the figure it is clear that 54, 2 is a better approximation to / than 
either 5s, i or Se,o- In fact, for this example, S$ t \ and 56, o lie close to each other and 
are each a “poor” approximation to /. We shall comment further on this example in 
section 8. 

If / is an even function of x, then results very similar to those above hold with each 
bk replaced by a^. In particular, for approximations of the form 

„ , X ^4o/2 + En=l cos(nx) 

SnM *> = i+e",c„ coster 

we find 

^ M 

Ak = O'k o ^ C m ((lk+m “f* Qk—m ) 1 0 — k — 

^ m = 1 




M — 771 

^ ] SkSm-\-k’> 1 < m , 


M 




9 J ’ 


0 


j = 0 


where the {s m } are now determined from the equations 


M 


(49) 


^ SmU]V+M+ 1—p— 77i — <2/V+M + l— p5 P 1? ”'7 




6. Analysis of the case M > 2. In this section, we prove some results analogous 
to those proved in section 4 for the case M = 1. 

THEOREM 4. Le/ / be an odd , 27 r-periodic, piece-wise smooth function and let its 
Fourier coefficients { 6 ^} satisfy the condition that = 0(1/ k p ), as k — ► oc, where p is a 
real , posziire number . Then , for any positive integer M > 1, C m = + 0(l/iV 2 ), /or 

1 < m < M, ir/iere c m> M are certain constants independent of A, and Bk = 0(1/A: P+2M ), 
as — > oc. More precisely, let b , ft, and 7 be constants (independent of k). Then if 

b k = (b/k p )( 1 + 0 /Jfc + 7 /Jk 2 + 0(1/ k 3 )), ask^oo, 


C m = 


Vm-W 


1 - 


m 2 (p + M 

2(2M 


i)(p + 2( M — 1)) 1 0(l/iV 3 ), 

f-l)JV 2 j V ' h 
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B k 


bpjp_+ 1 ) 


. ( T) -L A/f 


+ 0 ( 1/F +2M+1 ), 

as k,N — ► oo. Here y) = is the usual binomial coefficient. If f is an even, 

2-K-periodic, piece-wise smooth function, then these results hold with b k replaced by a k 
and B k replaced by A k . 

Proof: As in the proof of Theorems 1-3, we shall only outline the proof for the case 
when / is an odd, 27r-periodic, piece-wise smooth function. 

To see that the coefficients { C m } have the stated form as N — » oo, we first use the 
assumed form of the asymptotic form of the coefficients {6*} to write 


f501 bN+J = ( b /( N +j) P )( 1 +l 3 /( N + j) + 'l/( N +j) 2 + 0(l/N 3 )) 

= 6 c p (l +'(/? - jp)e + (7 -jfi(p+ 1) +j 2 p(p+ l)/2)e 2 + <3(e 3 )) , 

where we have defined e = 1/A. Next, we examine the behavior of the quantities {s m } 
as A — > 00. Using the linear equations (46), we solve for {s m } explicitly in terms of the 
coefficients {6*} for a few small values of M (see the Appendix). We then insert the 
expansions (50) into these explicit expressions for the {s m } and expand the resulting 
expressions for small values of e. In this way we obtain an expansion for each s m which 
is valid as c — + 0, i.e. as N — * 00. We then use induction to show that these expressions, 
for any positive integer M, can be written as 


Sm 


(- 1 ) 



p + M — 1 

N 


(51) M ( M ~ r ) il ± M ~ IHzM + ( M + 1 + - 2 ))p + rnf) 

V V 7 "- 1 / 2 pN 2 

for m = 1,2, ...,M, as N — *■ 00. Inserting equations (51) into (47) and expanding the 
resulting expressions for large values of A, we find that each coefficient C m has the 
asymptotic form indicated in the statement of the theorem. 

In a similar manner, using the expressions (50) with N replaced by k , along with 
the asymptotic form of the {C m }, in equations (48), we find that the coefficients {£?*} 
have the asymptotic form indicated in the statement of the theorem, as both k, N — ► 00. 
This completes the proof of our theorem. 

To illustrate these results, in Figure 7 we have plotted log \B k \j log k vs. 1 /logk 
for Example 1 with M = 1,2,3, and 4, for N = 20. The figure clearly illustrates how r 
much more rapidly the coefficients B k decay to zero (especially with increasing values 
of M) than the coefficients b k , which decay only like 1 /k. 

The results of Theorems 2 and 3 generalize to the case when M > 2 in a straightfor- 
ward manner. Using the results of Theorem 4, the proofs of the following two theorems 
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+ 0( l/Af 3 ), 



follow closely the proofs of Theorems 2 and 3, respectively, and, hence, we state them 
without proof. 

THEOREM 5. Let f be an odd, analytic, 2 k - periodic function and let its Fourier 
coefficients { 6 ^} satisfy the condition that bk — (b0 k fk p )( 1 + /? / k + 7 / k 2 + 0(1 / k 3 )) , as 
k — > 00 , where b, j3, and 7 are constants (independent of k), p is a real positive number, 
and 0 < |#| < 1. Then 

M — m Q2k f f M \ 

C m = 2(-1)"»“— /* M +0(1/I»), 

E“o (' ? ) 0“ 

and B k = O(\0\ k (l - 6 2 )(1 - k/N)/k p+M ) + O(\0\ k /N p + M+1 ), ask,N -*oo. Iff is an 
even, analytic, 2w-periodic function, then these results hold with bk replaced by ak and 
Bk replaced by A k - 

THEOREM 6 . Let f be an odd, 2i r-periodic, piece-wise smooth function and let 
its Fourier coefficients { 6 ^} satisfy the conditions of Theorem 1. Then the sequence of 
approximations converges to (f(x + ) + f(x~))j2, as N — ► 00 , for all —k < 

x < k. Moreover, the sequence {Sn,m(x)} converges like a Fourier series whose terms 
are 0(1/N P+2M ). These results also hold if f is an even, 2 k - periodic, piece-wise smooth 
function, if the Fourier coefficients {<Jjt} of f satisfy the conditions stated above with bk 
replaced by ak- 

If f is an odd, analytic, 2 k - periodic function and its Fourier coefficients satisfy the 
conditions of Theorem 2, then the sequence of approximations {SW.M^)} converges to 
f(x), as N — * 00 , for all —k<x<k. Moreover, the sequence converges like 

a Fourier series whose terms are 0(6 N / N p+M ). These results also hold if f is an even, 
analytic, 2 k - periodic function, if the Fourier coefficients {a*} of f satisfy the conditions 
stated above with bk replaced by ak- 

To illustrate Theorem 6 , in Figure 5 we have also plotted log \En,m\/ l°g |N| for 
Example 1 , with M = 2 and 3, as a function of l/log|iV|, for 3 < N < 50, with 
x = k/2. Here E n ,m(x) = (f(x + ) + f(x~))/2 - Sn,m(x). Using the ideas discussed after 
Theorem 3, the results presented in Figure 5 are consistent with the convergence rates 
of p = 5 (for M = 2) and p = 7 (M = 3), as predicted by Theorem 6 for the function 
considered in Example 1. 

7. Application. Applications of the results presented above to several classes of 
problems usually solved by Fourier series alone will be presented and discussed else- 
where. In this section, we present an application to a simple heat conduction problem 
to illustrate some of the potential of the method to improve the accuracy of approximate 
solutions obtained by partial sums of Fourier series. 

We consider the problem of determining the transient behavior of the temperature 
u(x,t) which satisfies the conditions 

u t = u xx , for 0 < x < k and t > 0 , 

u( 0 , t) — 0 = k( 7T , t ), for all t > 0 , and u(x, 0 ) = x/2, for 0 < x < k. 
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(52) 

with 


(In equation (52), the subscripts denote partial differentiation.) Using separation of 
variables, the solution is found in a straightforward manner to be 

(53) u(x ,1) = Jim where u<«»(x,i) = £ sin(ni) . 

This solution obviously has a discontinuity at x = n when t = 0. 

Using the formulas of section 2, we define a new class of approximate solutions 
{i/W 1 )} by 


(54) 


4) s EHLtS„(t)sin(nx) 
’ 1 + C\(i) cos(x) 


where the coefficients B n (t ) and C\{t) are defined by equations (13) and (14) with 6* 
replaced by ( — 1 ) k+1 e~ k2t /k. Thus, we find 


Ci(t) = 


2N(N + l)e“I 2 ^ + i) < 

(./V + 1 ) 2 + N 2 e~( 4N+2 ) t ’ 1-6 
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(55) 



Ci e~(* +1 ) 2< \ 1 

2 \ k- 1 + k + l ) J 


2 < k < N. 


In Figure 8 we have plotted u^.o) (dashed line) and u^ 3,1 ) (dotted line), along with the 
exact solution u(°°’ 0 ) (solid line) for t = 0.01 and t = 0.05. As the figure illustrates, id 3,1 ) 
is consistently a better approximation than u^ 4,0 l to the exact solution, and the quality 
of this approximation improves as t becomes larger. The corresponding comparison at 
t = 0 is the same as shown in Figure 1 for 0 < x < x. 


8. Conclusions and Discussion. We now make a few observations about the 
family of Fourier-Pade approximations discussed above, and also comment on several 
of their properties that need further investigation. 

First, we note that, using only the “information” contained in the first few Fourier 
coefficients of a function /, the functions Sn,m appear to provide new approximations 
to / which are consistently “better” than the Fourier partial sums Sjv, o °n which they 
are based. The sense in which these approximations are “better” can be interpreted in 
at least three different ways. First, Theorems 3 and 6 show that these approximations 
converge point-wise to (f(x + ) + f(x~))/2 at a faster rate than the original sequence 
{S N , 0 }, as the parameter N increases. Secondly, although the Gibbs phenomena is 
still present in the family {Sn,m} > the amplitude of the oscillations near a point of 
discontinuity of / appears to be mitigated, when compared with the oscillations present 
in Sff t Q. Finally, oscillations in 5 n,0i which lie outside a neighborhood of a point of 
discontinuity of / (or outside a region of large curvature of /), are noticeably damped 
in Sn Mi for M > 1, especially as N increases. These last two interpretations are 
illustrated in Examples 1 and 2. 

The fact that explicit expressions were derived for the coefficients A n , B n and C m , 
which appear in the definition of Sn,m, should be emphasized. (From a practical point 
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of view, it should be noted that it might be possible to develop an efficient, recursive 
algorithm, similar to algorithms presented by Baker [1], to compute these coefficients, 
since the equations which ultimately determine these quantities are recursive in nature. 
See, especially, equations (46)-(48), as well as the Appendix.) These explicit expressions 
not only facilitate the proofs of the stated theorems, but also illustrate some of the 
potential use of these approximations to construct approximate solutions to problems 
involving differential equations, especially partial differential equations. In particular, 
the Fourier coefficients {a n , 6 n } may be functions of one or more “other” variables, such 
as £, as in the example of section 7. The application of the basic ideas presented here to 
problems involving several different classes of partial differential equations is currently 
under investigation and will be reported elsewhere. 

As far as the mathematical properties of the family {•Sjv.Af} are concerned, only the 
most elementary properties have been investigated here. Many other questions, which 
are of both theoretical and practical importance, should be addressed. For example, for 
a fixed value of JV+M, which approximation S N , M is “best”? In Figure 9 we have plotted 
Se- q ,q, for q = 0 , 1 ,. ..,5 for Example 1 . In this case, the figure seems to suggest that 
perhaps 63,3 is “best” in some appropriate integral norm, although the improvement of 
the quality of the approximation of 83,3, say, over S^i is not dramatic. For this example, 
the function has essentially only one point of discontinuity (at x = tt) and hence Sn, 1 
appears to provide a “good” approximation. By contrast, the function considered in 
Example 3 has essentially two points of large curvature (near x = 8 and x = x - 6). 
For this example, 8 N ,i provided a poor approximation to /, while S N<2 yields a much 
improved approximation. More generally, the “best” Sjv.M will undoubtedly depend on 
both the general “shape” and smoothness properties of the function / , as well as on 
the particular norm used. 

From a practical point of view, a related issue concerns how best to represent a 
function with a discontinuity interior to [— x, x]. For example, consider the odd, piece- 
wise continuous, 2 x— periodic function / defined by f(x) = xx/ 2 , for — x /2 < x < x/ 2 , 
and f{x) = 0, otherwise, in [— x,x]. The Fourier coefficients {f> n } °f / are given by 
b n = (_i)(»-i)/2/ n 2, if n is odd, and b n = (x/(2n))(-l)<" +2)/2 , if n is even. The Fourier 
partial sum Sio,o associated with / is plotted in Figure 10. To approximate / by one of 
the functions Sjv.m, one possibility is to use the formulas above “blindly” and construct 
an approximation, say in a straightforward manner. This approximation is also 
plotted in Figure 10. As the figure clearly illustrates, neither Sio,i nor S 94 is a “good” 
approximation to /. However, we now observe that we can decompose S 2 n + 2,0 into the 
sum of two other partial sums, i.e. 


S 2 N+2,0 


^N+ 1,0 T ‘S'l'-f 1 , 0 ’ where 


SjiJJ .!, 0 = sin(*) - ^ sin(3x) + ^ sin(5x) - ..., 

•SjJi .,0 = l (lsin(2x) - jsin(4x) + jrsin(6x) - ...}. 
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We can now apply our formulas to each of these sums, separately, and construct two new 
approximations Sfy j and Sjgj. We then use 5^ j + as a new approximation to /. 

The approximation is also shown in Figure 10 and is clearly an improvement 

over either Sio^ or S^i . 

This example also serves to illustrate a possible advantage of the more general form 
of $n,m defined in equation (4). In particular, the only approximations that have been 
investigated here correspond to setting each D m = 0 in equation (4). More generally, the 
inclusion of both sine and cosine terms in the denominator of SW.m allows the possibility 
of “shifting” the location of an approximate “pole” of from either x = 0 or x = ±tt 

to an arbitrary point interior to the interval ( — tt, 7t). For example, when M — 1, the 
denominator of 5 jv,i can be expressed as 1 + Cj cos(x) + D x sin(x) = 1 + C x cos(x — <$). 
Thus, if Ci 1, for example, then the denominator will be small when x is near 6 
and, hence, Sn,i could potentially better simulate a function which has a singularity 
(or a “near” singularity) at x = 6. 

We also note that, while Theorems 3 and 6 establish the convergence of the sequence 
of approximations {5jv,Af} f° r a fixed value of M as N — * oo, the Gibbs phenomena 
has not been eliminated, although its effects seem to be mitigated. Thus, the rate of 
convergence of this sequence in regions near a point of discontinuity of / needs further 
study. For example, in the proof of Theorem 3, equation (35) establishes that, for any 
fixed value of x, with -7r < x < 7r, the terms r k are 0(l/k p+2 ) as k — > oo. However, the 
“practical” rate of convergence of the series is mitigated somewhat, especially near x = 
0, by the presence of the factor (1 — cos(x)) in the denominator of the expression for r k . 
This factor, in combination with the terms sin(Arx) in the numerator, also foreshadows 
the Gibbs phenomena that does remain. However, the fact that an explicit form for these 
terms is available should assist in the investigation of ways to improve convergence in 
these regions. This observation, as well as the other observations, questions, and several 
related issues raised here, are the subject of some current investigations and the results 
of these investigations will be reported elsewhere. 
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A. Appendix. In this appendix we present some explicit expressions for the quan- 
tities {s m } defined in section 5 in terms of the coefficients {6jt} by equations (46) and 
(49) when / is an odd function of x. Analogous expressions hold when / is an even 
function of x, with each bk replaced by a*. For M = 1,2, and 3, these quantities are 
given by 


M = 1 : «i = —bjv+i/bN] 


M = 2 : 


•si 


fo/v+lfyv ~ fr/V+2^N-l 

bN+ibN-i ~ b% 


, *2 = 


frjvfyv+2 — %+i . 
frjV+l^N-l - 


M = 3 : 


si = 


bN+l(b% ~~ f>/V- lV+l) + f>N+ 2 (^Y+l^Y -2 ~ f>N-l^jv) + f>N+ 3(^-1 ~ fyvfyv- 2) 


det 


«2 - 


6/V+l (6/v+2^/V-l — ~H fojV+2(frjy ~ ^N+2^N-2) ~t~ ^N+3(^N+l^iV-2 bfi/biN- -l) 


det 


«3 


6 Af +1(^+1 — 6 ^^+ 2) + ^+2(^+2^-! ~ bffbpf+i) + f>N+ 3 (^ ~ ^N+l^lV-l) 


det 


det = 6jv_ 2( — + b^b^+2) + b^-i(—bN-ibN+2 + 2&jv6/v+i) + ^v( — 
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Figure 1: A plot of the function f(x ) (solid line), the approximations 5 4)0 (x) 
(dashed line) and 5*3,1 (x) (dotted line) for Example 1. Note that the oscillations present 
in 5 4 , o away from the singularities of / have been virtually eliminated in 5 31 . 
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Figure 4: A plot of jfc 2 |A*|/|a*| (circles) for Example 2 as a function of 1/k, for 
2 < k < N, with 3 < N < 20. The figure illustrates that the coefficients {A k } decay 
faster than the coefficients {a^}, asymptotically by a factor of k 2 , as predicted by 
Theorem 2. As in Example 1, an interesting asymptotic behavior of A N - q , as N — ► oo 
for a fixed value of q , is also predicted. 





Figure 5: A plot of log \E NM \f log \N\ for Example 1 as a function of 1 / log |7V| for 
3 < N < 50, at x = 7r/2, for M = 0, 1,2, and 3. Here E NM (x) = (f(x+) + f(x~))/2 - 
Sn,m{x)- We note that an intercept of 2 M + 1, as predicted by Theorems 3 and 6, is 
clearly consistent with these plots. 






Figure 6: A plot of the function f(x) (solid line), the approximations S 6 ,o(*) 
(longer dashed line), 5^1 (x) (shorter dashed line), and S^i^) (dotted line) for Example 
3 . Note that the oscillations present in Se,o away from the regions of large curvature 
of / have essentially been eliminated in S^t while is not much of an improvement 
over S e ,o- 
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Figure 7: A plot of log|5^|/logA: vs. l/log& for Example 1 with M = 1,2,3, 
and 4, for JV = 20. The increased rate of decay to zero of the coefficients with 
increasing values of M is clearly evident, especially when compared with original Fourier 
coefficients 6*, which decay only like 1 jk. 





Figure 8: A plot of the exact solution u (oo,0) (solid line), as well as the approxima- 
tions u (4 ’ 0) (dashed line) and id 3 - 1 ) (dotted line), for the simple heat conduction problem 
of section 7, when t = 0.01 (upper set of curves) and t = 0.05 (lower set of curves). The 
corresponding comparison at t = 0 is the same as shown in Figure 1 for 0 < x < 7r. 
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e approximations i> 6 - 9)9 , ior 
;ems to suggest that perhap 
= 6 . 





Figure 10: For the 27 r-periodic function / defined by f(x) = kx/2, for — 7 t/2 < 
x < 7r/2, and /(x) = 0, otherwise, in [— 7r,7r], a plot of f(x) (solid lines) and the 
approximations 5io,o (longer dashed line), Sg.i (shorter dashed line), and + S[ 2 J 
(dotted line) for 0 < x < ic. The improved quality of the approximation + S{ 2 }, 
over either Sio i or Sg i, is apparent. 
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